Oscillatory Universes in Loop Quantum Cosmology and Initial Conditions for Inflation 
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Positively-curved, oscillatory universes are studied within the context of Loop Quantum Cosmol- 
ogy subject to a consistent semi-classical treatment. The semi-classical effects are reformulated in 
terms of an effective phantom fluid with a variable equation of state. In cosmologies sourced by a 
massless scalar field, these effects lead to a universe that undergoes ever-repeating cycles of expan- 
sion and contraction. The presence of a self-interaction potential for the field breaks the symmetry 
of the cycles and can enable the oscillations to establish the initial conditions for successful slow-roll 
inflation, even when the field is initially at the minimum of its potential with a small kinetic energy. 
The displacement of the field from its minimum is enhanced for lower and more natural values of 
the parameter that sets the effective quantum gravity scale. For sufficiently small values of this 
parameter, the universe can enter a stage of eternal self-reproduction. 



I. INTRODUCTION 

Cyclic/oscillatory universes have a long history in cos- 
mology 1]. Originally, one of their main attractions 
was that initial conditions could in principle be avoided. 
Closer consideration of such models, however, revealed 
severe difficulties in their construction within the con- 
text of general relativity (GR). Apart from entropy con- 
straints, which restrict the number of bounces in the 
past, the central difficulty is that when treated classi- 
cally any bounce would be singular, thereby resulting 
in the breakdown of GR. Recent developments in M- 
theory inspired braneworld models have renewed inter- 
est in cyclic/oscillatory universes 0, although problems 
still remain in such models in developing a successful 
treatment of the bounce. An oscillating universe that ul- 
timately undergoes inflationary expansion after a finite 
number of cycles has also been investigated Q. However, 
a physical mechanism for inducing the bounces was not 
employed in this model. 

Our aim here is to study oscillatory universes within 
the context of Loop Quantum Cosmology (LQC) which 
is the application of Loop Quantum Gravity (LQG) to 
an homogeneous minisuperspace environment. LQG is 
at present the main background independent and non- 
perturbative candidate for a quantum theory of gravity 
(see for example [j, |5| ) • This approach provides a (dis- 
crete) description of high energy dynamics in the form 
of a difference equation. An important consequence of 
this discretization is the removal of the initial singularity 
J6J. As the universe expands and its volume increases, it 
enters an intermediate semi-classical phase in which the 
evolution equations take a continuous form but include 
modifications due to non-perturbative quantization ef- 
fects 0- It has recently been shown that a collapsing 
positively-curved Friedmann-Robertson- Walker (FRW) 
universe containing a massive scalar field can undergo 
a non-singular bounce within the semi-classical region 
8J. Other studies have shown that for a flat universe 
semi-classical LQC effects increase the parameter space 
of initial conditions for successful inflation, and this be- 



havior is found to be robust to ambiguities in the theory 

MM 

We study the behaviour of a positively-curved universe 
sourced by a massless scalar field and show that the com- 
bined effects of the LQC corrections and the curvature 
enable the universe to undergo ever-repeating cycles. The 
presence of a self-interaction potential for the field breaks 
the symmetry of the cycles and can establish the condi- 
tions for inflation, even if the field is initially located in 
the minimum of its potential. The viability of this mech- 
anism is studied by identifying the constraints that need 
to be satisfied for a consistent semi-classical treatment. 



II. EFFECTIVE FIELD EQUATIONS IN LOOP 
QUANTUM COSMOLOGY 

We consider positively-curved FRW cosmologies 
sourced by a scalar field <f> with self-interaction poten- 
tial V(4>). The semi-classical phase of LQC arises when 
the scale factor lies in the range at < a < a*, where 
a,i = y/jlpi, a* = y/jj/Slpi, 7 = ln2/\/37r « 0.13 and j 
is a quantization parameter which must take half integer 
values. Below the scale Oj, the discrete nature of space- 
time is important, whereas the standard classical cos- 
mology is recovered above a* . The parameter j therefore 
sets the effective quantum gravity scale. The modified 
Friedmann equation is given by 
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where the quantum correction factor D(q) is defined by 

D(q) = QLy^/2{ 7 [(g + l)"/*-|«-l| U/4 " 

"(q + I) 7 / 4 - sgn(g - l)|g - 1| 7 / 4 ] } 6 (2) 
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with q = (a/ a*) 2 . Eq. J2J) represents an approximate ex- 
pression for the eigenvalues of the inverse volume opera- 
tor 0| ■ As the universe evolves through the semi-classical 



phase, this function varies asflot a 15 for a <C a*, has 
a global maximum at a ~ a* and falls monotonically to 
D = 1 for a > a*. 

The scalar field equation has the form 
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where a prime denotes differentiation with respect to the 
scalar field. Differentiating Eq. QJ and substituting Eq. 
P)l then implies that 
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Equations JQ and can be written in the standard 
form of the Einstein field equations sourced by a perfect 
fluid: 
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Pcff = -3i/(/9 cff +p e ff) , 
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define the effective energy density and pressure of the 
fluid, respectively. The effective equation of state, w = 
Pes/Pefi, is given by 
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In this picture the effects of the LQC corrections on 
the cosmic dynamics are parametrized entirely in terms 
of the equation of state. When d In D/d In a > 6, the fluid 
represents 'phantom' matter (w < —1) that violates the 
null energy condition (p c g + p c g > 0) and, numerically, 
this is equivalent to a < a p h = 0.914a*. This condition 
depends only on the quantization parameter j and is in- 
dependent of the potential. Since, within the context of 
GR, a violation of the null energy condition can lead to 
a non-singular bounce, this provides an alternative de- 
scription of how LQC effects can result in a bouncing 
cosmology when the universe contracts below a critical 
size: the decreasing energy density of the phantom fluid 
during the collapsing phase is eventually balanced by the 
growing curvature term in the Friedmann equation. 

In the following Sections, we develop this picture fur- 
ther for the case of massless and self-interacting scalar 
fields, respectively. 



III. OSCILLATIONS WITH A MASSLESS 
SCALAR FIELD 

The field equation J3J) for a massless scalar field (V = 
0) admits the first integral: 
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where D; n it = ^?(ainit/o*) an d a subscript 'init' denotes 
initial values. Without loss of generality we choose initial 
conditions such that a» < a; n it < a p h, </>i n it > and 
-ffinit = 0. The subsequent cosmic dynamics can then be 
divided into four phases (Fig. QJ. 

Phase I: The universe is effectively sourced by phantom 
matter since a < a p h- This drives an epoch of superinfla- 
tionary expansion (H > 0) until din D/d In a = 6. The 
scalar field accelerates during this phase. 

Phase II: The expansion rate now slows down. The 
scalar field continues to accelerate until the logarithmic 
slope of the quantum correction function has fallen to 
din D/d In a = 3. For a > a*, the energy density of the 
field, Pff, ex 1/a 6 , falls more rapidly than the curvature 
and the expansion eventually reaches a turnaround. 

Phase III: Eq. @ implies that the scalar field begins 
to accelerate immediately after the turnaround and con- 
tinues to do so until the universe has collapsed to the 
point where d In D/d In a = 3 once more. Shortly after- 
wards, the universe contracts below the critical scale o p h, 
at which point the Hubble parameter begins to increase. 

Phase IV: Substituting Eq. (UHJ into Eq. Q im- 
plies that the effective energy density varies as p e g = 
<j> 2 /(2D) oc D/a 6 . It is therefore decreasing during this 
phase, since dln£)/dlna > 6, whereas the curvature is 
growing. Hence, a scale is reached where the Hubble 
parameter vanishes instantaneously. Eqs. QJ and Q 
imply that the second time-derivative of the scale factor 
is positive during this phase and the universe therefore 
undergoes a non-singular bounce into a new phase of su- 
perinflationary expansion (Phase I). 

The effective equation of state (0 is independent of 
the field's kinetic energy when V = 0. This implies that 
identical cycles are repeated indefinitely into the future 
(and the past), as summarized in Fig.[2Ji. During the col- 
lapsing phases, the field retraces the trajectory it mapped 
out during the expanding phases. Oscillatory behaviour 
is therefore possible in the absence of an interaction po- 
tential. Furthermore, Eq. (|10fl implies that the field's 
kinetic energy never vanishes during the cycle since the 
scale factor remains finite. The value of the field there- 
fore increases monotonically with time. 



IV. SELF-INTERACTING SCALAR FIELD 

We now consider how self-interactions of the scalar 
field modify this cyclic dynamics. We make very weak 
assumptions about the potential, specifying only that it 
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FIG. 1: Time evolution of the scalar field <j>, the logarithmic 
scale factor and the Hubble parameter when the potential 
V — 0, j — 100, ainit/a* = 0.9 and Hi = 0. Axes are labeled 
in Planck units. 
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FIG. 2: (a) Schematically illustrating the logarithmic varia- 
tion of the effective energy density of a massless scalar field 
(solid line) and the curvature term in the Friedmann equation 
(dashed line) . The universe oscillates indefinitely between the 
intersection points of the two lines, (b) Schematically illus- 
trating the effects of introducing a self-interaction potential 
for the field. The cycles are eventually broken as the potential 
becomes dynamically significant, thereby resulting in slow- 
roll inflation. In both figures, the slope of the trajectories is 
given by dlnp c ff/dlna = — 3[1 + w(a)]. 



has a global minimum at VJ T1 i n (0) = and is a positive- 
definite and monotonically varying function when (f> =/= 
such that V" > 0. We suppose the field is initially located 
at the minimum of its potential with the same initial 
conditions as those considered in Section fff . 

The qualitative dynamics of the universe is illustrated 
in Fig. |2h- In general, the path of the field in the 



{lnpeff, lna} plane is determined by the variation of the 
effective equation of state ©. The gradient of the trajec- 
tory is given by d\np e g/d\na = — 3[l + w(a)] and there is 
a turning point whenever the scale factor passes through 
a p h, regardless of the form of the potential. 

Two factors determine how the equation of state 
changes during each cycle and these can be parametrized 
by defining the quantities B = [(6 — dlnD/dhia)/6] and 
C = 20 2 /[0 2 + 2DV], respectively. For a massless field, 
C = 2, and the deviation of the equation of state away 
from the value w = — 1 is determined by the magnitude 
of B. The effects of the potential are parametrized by C. 
Introducing a potential energy into the system necessar- 
ily implies C < 2 and it follows that for a given value of 
the scale factor, the equation of state is closer to —1 than 
in the massless case. 

Since the field moves monotonically up the potential, 
the gradient term DV' in the scalar field equation © 
is more significant in a given cycle relative to the previ- 
ous one. This implies that the scalar field gains kinetic 
energy less rapidly during Phases I and III and loses ki- 
netic energy more rapidly during Phases II and IV. Con- 
sequently, the value of C at the end of Phase I of a given 
cycle is smaller than the value it had at the corresponding 
point of the previous cycle. In other words, the rate of 
change of the trajectory's gradient around a p h becomes 
progressively smaller with each successive cycle. 

A massless field begins accelerating immediately after 
the turnaround and bounce have been attained. How- 
ever, Eq. © implies that <\> — —DV' < at the in- 
stant when the Hubble parameter vanishes, so the field 
does not accelerate immediately after the turnaround 
(bounce) ; there is a short delay during which C continues 
to decrease. The trajectory of the field immediately be- 
fore the turnaround (bounce) is therefore slightly steeper 
than the trajectory immediately after, i.e., the trajectory 
for Phase III (I) lies below that of Phase II (IV). More- 
over, because the field does not acquire as much kinetic 
energy during Phase III as it lost during Phase II, the 
value of C (for a given value of the scale factor) is smaller 
during Phase III than it was during Phase II. Thus, the 
Phase III part of the trajectory always lies below that 
of Phase II and, similarly, the trajectory of Phase I lies 
below that of the Phase IV trajectory of the previous cy- 
cle. The potential therefore breaks the symmetric cycles 
of the massless field and the effective energy density at 
the end of Phases I and III falls with successive cycles. 

We have implicitly assumed that the kinetic energy of 
the field never vanishes during Phase II. However, the po- 
tential becomes progressively more important with each 
completed cycle and this assumption must necessarily 
break down after a finite number of cycles have been 
completed. Since the effects of the quantum correction 
function are negligible once the universe has expanded 
beyond a*, the energy density of the field during Phase 
II redshifts more rapidly than the curvature term if the 
strong energy condition is satisfied, i.e., if (j) 2 > V. As 
the field moves monotonically up the potential, it be- 



comes progressively harder to maintain this condition. 
Eventually, therefore, a cycle is reached where this con- 
dition is violated during Phase II. This leads inevitably to 
an epoch of slow-roll inflationary expansion, as the field 
slows down, reaches a point of maximum displacement 
and moves back down the potential [jj|. The transition 
into a potential-driven, slow-roll inflationary epoch is 
shown in Figs. |2b and [3] 

Thus far, we have discussed the cyclic dynamics for 
a given value of the quantization parameter, j. In spa- 
tially flat models, the maximum value attained by the 
field as it moves up the potential increases for higher val- 
ues of j, since the universe undergoes greater expansion 
before the end of the supcrinflationary, semi-classical 
phase. For positively-curved models, on the other hand, 
the field moves further up the potential before the onset 
of slow-roll inflation for smaller values of j. This be- 
haviour can be understood qualitatively by making some 
simplifying assumptions. Let us assume that the transi- 
tion from Phase I to II occurs at a — a»; that D ex a 15 
for a < a*; that D = 1 for a > a*; and that the poten- 
tial becomes dynamically significant only during the last 
cycle before slow-roll inflation. At this level of approxi- 
mation, the field's kinetic energy varies as 2 |a<a, oc a 24 



and 4> \ a >a, oc a , respectively, so • 
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transition. At the turnaround (a = a t ), </> 2 = 2 (a*/at) 6 , 
but since <j) 2 rs 1/a 2 , it follows that </> 2 oc j~ 15 ^ 2 . Hence, 
the kinetic energy of the field at the turnaround is higher 
for lower values of j and, consequently, slow-roll infla- 
tion occurs (V > 4> 2 ) when the field is further from the 
minimum of the potential. 

To summarize, the cyclic nature of positively-curved 
cosmologies within LQC can in principle establish the 
conditions for slow-roll inflation even if the inflaton is 
initially located in the minimum of its potential. We 
now proceed to consider the parameter space of a viable 
model of inflation. 
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FIG. 3: Illustrating the time evolution of the scalar field cj> 
(top panel) and the logarithmic scale factor (middle panel) 
for a quadratic potential V = m 2 (f> 2 /2. The bottom panel 
illustrates how the effective kinetic energy <j> /2D (solid line) 
and potential energy (dashed line) of the field vary compared 
to the curvature term (grey line) . The cycles end and slow-roll 
inflation commences at the point when the potential begins 
to dominate the kinetic energy. Initial conditions are chosen 
such that </>i n it = -ffinit = and ainit/a* = 0.9, with m = 
10 — 6 Zp 1 1 and j — 5 x 10 11 (see the text for details). The axes 
are labeled in Planck units. 



Friedmann equation IjTjl: 
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On the other hand, constraint (b) results in an upper 
bound on the field's kinetic energy. Eq. Q implies that 
\H\ is maximized when dlnD/dliia w 6 and substituting 
Eq. i|10|) into the Friedmann equation JJ| then implies 
that |-ffmax|ai < 1 for 



V. PARAMETER SPACE OF A VIABLE MODEL 

A number of constraints must be satisfied by any suc- 
cessful inflationary scenario of the type outlined above. 
In particular, the field must be sufficiently displaced from 
the minimum of its potential at the end of the oscilla- 
tory phase for the horizon problem to be solved. Typ- 
ically this requires at least 60 e-folds of accelerated ex- 
pansion. Self-consistency of the semi-classical analysis 
also requires that: (a) H 2 must be non-negative; (b) the 
Hubble length must be larger than the limiting value of 
the scale factor, \H |oi < 1; and (c) the scale factor at the 
bounce must exceed ai . Imposing constraint (b) is equiv- 
alent to requiring that energy scales during the classical 
regime do not exceed the Planck scale J3- 

Constraint (a) necessarily implies that the initial ki- 
netic energy of the field is bounded from below by the 
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where Z? p h = D(a p h/a*)- The discussion leading to Fig. 
|2b has shown that the maximum value of the square of 
the Hubble parameter during a given cycle is smaller than 
that of the previous cycle. Thus, if condition l|12fl is sat- 
isfied in the first cycle, where the potential is negligible, 
|.ff|aj < 1 during all subsequent cycles. 

Fig. 01 illustrates the region of parameter space where 
constraints (|llfl and l|12|) are simultaneously satisfied. 
These constraints are independent of the potential. In- 
creasing the value of j reduces the upper and lower 
bounds on the initial kinetic energy but widens the range 
of allowed values of ainit/a*- A further consequence of 
Fig. 21 is that for a given value of j there exists an upper 
limit to the size of the universe at the turnaround that 
is consistent with the semi-classical dynamics. In other 




FIG. 4: Constraints on the initial velocity of the field |</>|mit 
for a given initial value of the scale factor, ttinit/a*, for dif- 
ferent values of the quantization parameter, j. The axes are 
labeled in Planck units and the shaded areas represent the 
regions where constraints Hill and 111 21 are satisfied. For all 
values of j, the area of the shaded region is finite, and the 
points of intersection occur further from amit/a* ~ 1 as j is 
increased. Note that for amit > a*, the universe is initially 
in a contracting phase and subsequently bounces, after which 
the behaviour discussed in the text is followed. 



words, if the universe is too large at the turnaround, the 
kinetic energy of the scalar field will exceed the Planck 
scale before the universe has contracted below a*. This 
implies, in particular, that if the inflaton did not decay 
at the end of inflation, future cycles after the slow-roll 
inflationary epoch could not be realized within this semi- 
classical framework. 

Constraint (c) is not necessarily satisfied if the bounce 
occurs at progressively smaller values of the scale factor. 
It is necessary, therefore, to consider this constraint for 
each specific model. We have considered a quadratic po- 
tential V(4>) = m 2 (j) 2 /2, where the mass of the field is 
set by the COBE normalization to be m = 10 — 6 Zp 1 1 . We 
have verified numerically that constraints (a)-(c) remain 
satisfied for all initial values of the parameters contained 
within the shaded regions of Fig. 0] and, furthermore, 
that the maximum value attained by the field is suffi- 



cient to solve the horizon problem, i.e., </> max > 3Zpj . 
Indeed, for the initial condition <Zi n it = 0.9a* considered 
in Fig. |21 sufficient inflation is possible for j < 5 x 10 11 . 
Numerical integration indicates that successful slow-roll 
inflation is also possible for other potentials such as quar- 
tic and hyperbolic models. 



VI. CONCLUSION 

We have shown that loop quantum gravity corrections 
to the Friedmann equation in a positively-curved uni- 
verse enable a scalar field to move up its potential due 
to a series of contracting and expanding phases in the 
cosmic dynamics. The oscillatory nature of such mod- 
els can be described by reformulating the semi-classical 
dynamics in terms of an effective phantom fluid. In gen- 
eral, the field is able to move further from the potential 
minimum for lower and more natural values of the quanti- 
zation ambiguity parameter j. This is important because 
for sufficiently small j, the field may move far enough up 
the potential before the cycles are broken for the universe 
to enter a stage of eternal self-reproduction [llj . Eternal 
inflation occurs if 3 V 2 < 1287WpiV 3 and for a quadratic 
potential this implies that > l/2m 1 / 2 /|, 1 w 500/pj 1 [i"l| . 
Numerical simulations indicate, for example, that this 
condition is satisfied for j < 2 x 10 7 when ai n ; t = 0.9a*. 
We conclude, therefore, that even if the field is at the 
minimum of its potential, a wide range of initial condi- 
tions leads to successful (and eternal) inflation. Since the 
assumptions made about the form of the potential were 
weak, this mechanism should be very generic. 

Finally, the oscillatory dynamics also arises when the 
minimum of the potential is at V min < 0. In positively- 
curved universes, the same mechanism enables the field 
to work its way out of the negative region and continue 
up the potential until the cycles are broken, as discussed 
in Section IV. For spatially flat universes, the picture is 
similar although the oscillations come to end as soon as 
the field reaches a positive region of the potential. This is 
particularly interesting given that negative potentials are 
known to arise in string/M-theory compactifications (e.g. 
|P2j|). It would be interesting to explore these possibilities 
further. 
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